Interference in the resonance fluorescence of two incoherently coupled transitions 
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The fluorescence light emitted by a 4- level system in J = 1/2 to J = 1/2 configuration driven by 
a monochromatic laser field and in an external magnetic field is studied. We show that the spectrum 
of resonance fluorescence emitted on the n transitions shows a signature of spontaneously generated 
interference effects. The degree of interference in the fluorescence spectrum can be controlled by 
means of the external magnetic field, provided that the Lande g-factors of the excited and the 
ground state doublet are different. For a suitably chosen magnetic field strength, the relative weight 
of the Rayleigh line can be completely suppressed, even for low intensities of the coherent driving 
field. The incoherent fluorescence spectrum emitted on the n transitions exhibits a very narrow peak 
whose width and weight depends on the magnetic field strength. We demonstrate that the spectrum 
of resonance fluorescence emitted on the a transitions show an indirect signature of interference. A 
measurement of the relative peak heights in the spectrum from the a transitions allows to determine 
the branching ratio of the spontaneous decay of each excited state into the a channel. 

PACS numbers: 42.50.Ct, 32.50.+d, 42.50.Lc, 42.50.Xa 



I. INTRODUCTION 

Since the emergence of quantum mechanics, quantum 
interference has been regarded as one of the most exciting 
and intriguing aspects of quantum theory . Although 
interference effects are present in almost all areas of quan- 
tum mechanics, some of them particularly attracted the 
attention of many scientists. In the following, we will give 
two examples of physical systems that are well known in 
the context of interference effects and are both related to 
the work presented here. 

First of all, we would like to mention the so-called V- 
system that has been intensively discussed by theoret- 
ical means. This atomic level scheme is comprised of 
two near-degenerate excited levels and one ground state, 
and many authors demonstrated that a rich variety of 
interference effects should be observable in this system. 
These effects include the modification and quenching of 
spontaneous emission [3- 0- [J ■ and several schemes to 
control spontaneous emission by means of external fields 
have been suggested |H SHI- Furthermore, it has 
been shown that quantum interference leads to strong 
modifications of the spectrum of resonance fluorescence, 
and for suitable parameters the complete supp ression 
of resonance fluorescence is achievable 0, 0, • 
The emitted fluorescence light also displays highly non- 
classical features like extre mely strong intensity-intensity 
correlations and squeezing |14L |X5| . 

However, all these schemes rest on the existence of 
spontaneously generated coherences between the two up- 
per levels that can only arise if the dipole moments be- 
tween the two upper and the lower level are parallel or 
at least non-orthogonal. This requirement is very hard 
to meet in an experiment, since appropriate atomic sys- 
tems are not known up to now. In order to circumvent 
this problem, an exp eriment with a molecular system has 
been performed but the experimental results could 
not be reproduced yet 01 • A recent experiment demon- 
strates the existence of spontaneously generated coher- 



ences between spin states in quantum dots |Tsj . 

One of the most famous interference effects is certainly 
Young's double-slit experiment, especially because it al- 
lows to explore fundamental concepts of quantum me- 
chanics such as the principle of complementarity in a 
very simple setup. Celebrated thought experiments like 
Feynman's light microscope |19| and Einstein's recoiling 
slits |2(j| | employ the position-momentum uncertainty re- 
lation to demonstrate that it is impossible to observe the 
wave and the particle nature of the interfering quantities 
(for example, electrons or photons) at the same time. In 
recent years, a prop osal b y Scully et. al. [2l| gave rise to 
a lively debate Hi |H Hf |H M, IH E| on the interre- 
lation between the principle of complementarity and the 
position-momentum uncertainty relation. 

A beautiful realization of Young's two-slit experiment 
was performed by Eichmann et. al. 29] and subsequently 
discussed by several authors [13, l3ll l3^ | . In this ex- 
periment, the slits are represented by two 198 Hg + ions 
in a trap that are irradiated by a coherent laser field, 
and the interference pattern formed by the scattered 
light was observed. The level scheme of each of the 
two atoms can be modeled by a J = 1/2 to J = 1/2 
transition that is also in the focus of the work presented 
here; a schematic representation of this four-level sys- 
tem is shown in Fig. ^ 33], 02 HH- The transitions 
1 1 ) *-> | 4 ) and | 2 ) <-> | 3 ) couple to er+ and <j~ polarized 
light, respectively, and will be referred to as the a tran- 
sitions. By contrast, the tt transitions | 1) <-> |3) and 
| 2 ) «-> | 4 ) couple to light linearly polarized along e z , 
and their dipole moments are anti-parallel. The four-level 
system of Fig. ^ is thus a realistic level scheme with non- 
orthogonal dipole moments which can be found in real 
atoms. However, it cannot be expected that this four- 
level system displays the same interference effects that 
were predicted for the V-system with parallel dipole mo- 
ments since there is a striking difference between them. 
In the case of the V-system, both transitions from the 
upper levels end up in the same ground state, while the 
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FIG. 1: Schematic representation of the four-level atom of in- 
terest. The two upper and lower levels are Zeeman sub-levels 
with m 3 = ±|. Each upper state can decay by a dipole al- 
lowed transition to both ground states. The Zeeman splitting 
of the magnetic sub-levels is not to scale. 



two 7T transitions of our four-level system start and end 
up in different states that are orthogonal to each other. 
Thus the question arises whether interference effects can 
also be observed in single J — 1/2 to J = 1/2 systems. 

Recently we investigated the fluorescence light emitted 
on the 7r transitions in the case of the degenerate system 
(B = 8 = in Fig. ^| and for a monochromatic driving 
field polarized along e z . It has been shown (3(| that the 
spectrum of resonance fluorescence indeed exhibits a sig- 
nature of vacuum-mediated interference effects, whereas 
the total intensity is not affected by interference. It has 
been demonstrated that this result is a consequence of the 
principle of complementarity, applied to time and energy. 

Here we generalize our analysis to the non-degenerate 
system depicted in Fig. ^ an d include the fluorescence 
light emitted on the a transitions in our discussion. Since 
it is possible to discriminate between the fluorescence 
light that stems from the n- and the a transitions simply 
by means of a polarization-dependent detection scheme 
(Seem, these two contributions will be discussed inde- 
pendently in Sees. Hill and IIVI respectively. 

We find that for the spectrum from the 7r transitions, 
the degree of interference in the coherent and incoherent 
part of the spectrum strongly depends on the frequency 
difference 6 between the two tt transitions. The spectral 
properties of the fluorescence light can thus be controlled 
by means of an external magnetic field that determines 
the value of 5. For example, the relative weight of the 
Rayleigh line can be completely suppressed for a suit- 
able value of <5, even for low intensities of the coherent 
driving field (Sec. IIII A*) l. The dependence of the incoher- 
ent spectrum of resonance fluorescence on the parameter 
8 is discussed in Sec. IIII Bl Section IIII CI demonstrates 
how the interference terms alter the fluorescence spec- 
trum emitted on the 7r transitions for different regimes of 
the driving field strength. The experimental observation 
of the spectra including the interference terms could pro- 



vide evidence for vacuum-mediated interference effects in 
an atomic system. 

The fluorescence spectrum emitted on the a transi- 
tions of the degenerate system is discussed in Sec. IIVI It 
only consists of an incoherent part and shows an indirect 
signature of interference since the relative peak heights 
depend on the interference terms. The measurement of 
the relative peak heights would also allow to determine 
the branching ratio of the spontaneous decay of each ex- 
cited state into the a channel. 

Section IV1 provides a detailed discussion of our results. 
In our our previous work |36j , we interpreted the interfer- 
ence effect in the spectrum of resonance fluorescence in 
terms of interferences between transition amplitudes that 
correspond to different time orders of photon emissions. 
Here we support this explanation by a formal argument, 
and the continuous transition from perfect frequency res- 
olution to perfect temporal resolution is studied in more 
detail. For a weak driving field and certain values of 
the parameter 5, the incoherent spectra of the n and a 
transitions contain a very narrow peak whose width is 
smaller than the natural linewidth. We explain these 
narrow structures in terms of electron shelving. Finally, 
a brief summary of our results is given in Sec. IVII 



II. EQUATION OF MOTION AND DETECTION 
SCHEME 



We now return to the level scheme in Fig. ^ Note 
that we allow the Zeeman splitting of the excited and 
the ground state magnetic sub- levels to be different, since 
the Lande g-factors will not necessarily be the same for 
these two multiplets. For example, in the case of the 
6s 2 5i/2 — Qp 2 Pi/2 transition in 198 Hg + the g-factor for 
the excited states is given by 2/3, and for the ground 
states it takes on its maximum value of 2. The matrix 
elements of the electric-dipole moment operator d can be 
found from the the Wigner-Eckart theorem [37| and are 
given by 

d x = < 1 3) = —^=Ve z , d 2 = (2|d|4) = -d! 

d 3 = (2|d|3) = y|p e (-), d i = (l\d\4)=d;. (1) 

In this equation, the circular polarization vector is de- 
fined as e' - ) = (e x — i e y ) and T> denotes the re- 
duced dipole matrix element. We assign to each of the 
four dipole-allowed transitions a resonance frequency uJi 
(i € {1,2,3,4}). If the splitting between the magnetic 
sub- levels vanishes (i.e. B = 8 = 0), these four frequen- 
cies are equal. 

We are interested in the time evolution of our four level 
system driven by a monochromatic field of frequency 
that is linearly polarized along the z axis, 



E(t) = E Q e 



-lU)L,t 



(2) 
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and c.c. stands for the complex conjugate. With this 
choice of polarization, the electric field couples only to 
the two anti-parallel dipole moments d\ and d.2- In the 
rotating wave approximation, the interaction Hamilto- 
nian takes the form 

V = (ilis- A 2i )hVLe- iuJLt + h.c. , (3) 

where the atomic transition operators are defined as 
■Aij = I * ) ( 1 1 an d the Rabi frequency is given by 
Q = EqV I \\fZK). The atomic Hamiltonian can be writ- 
ten as 

H = hujx A X1 + h{uj 2 + B) A 22 + hB A44 , (4) 

where w\ stands for the resonance frequency of the 1 <-> 3 
transition and u>2 = wi + S is the resonance frequency on 
the 2 <-> 4 transition. In a rotating frame defined by the 
unitary transformation 

W = exp [( A n + A 22 ) iu L t] , (5) 

the master equation for the density operator g = WgW^ 
reads 

1 

Q = -r [H,g] + C 1 g . (6) 

In this equation, the Hamiltonian is given by 

H = -h[AA n + {A-5)A22-B(A22+A 44 )] + 

[{Ax3-A M )m + h.c.], (7) 

A = u>i — u>% is the detuning of the driving field from 
resonance with the 1 «-> 3 transition, A— 15 is the detuning 
on the 2 «-> 4 transition and the damping term £ 7 g takes 
the form 

I 2 

c iQ = E 7 'y { i s i~> s j $ + [e s t> s j ] } 

-yE + ■ (8) 

i=3 

The transition operators 5^ are defined as 

5^ = A13 , s£ = A24 , = A23 , S4" = A14 , (9) 

and 5~ = (S^y. The decay constant on each of the 
cr transitions is denoted by j a , the parameters 7y are 
determined by 

^ = WW\^ *.ie{i,2}, (10) 

and 71 and 72 are the decay constants of the 7r tran- 
sitions (see Fig. P). For i = j, the parameters 7^ are 
equal to the decay rates of the 7r transitions, 711 = 71 
and 722 = 72- Although 71 and 72 are equal in our 



setup, we will continue to label them differently to facil- 
itate the physical interpretation later on. Since d\ and 

d, 2 are anti-parallel, the cross-damping terms are given by 
712 = 721 = — v/7i72- These terms allow for the possibil- 
ity of coherence transfer from the excited to the ground 
state doublet. 

The decay rates 71, 72, 7 CT can be related to the total 
decay rate 7 = 71 +"f a — 72 +7o- of each of the two excited 
states through the branching probabilities and b a , 

71 = 72 = ^7 and 7c = &<r7 ■ (H) 

According to the Clebsch-Gordan coefficients, we have 
&7r = 1/3 and b a — 2/3. Although we will keep the 
symbols b^ and b a in formulas, we will always assume 
these values whenever a concrete evaluation is performed, 

e. g. in figures. 

Next we employ the normalization condition Tr(g) = 1 
to eliminate the matrix element 044 from the master 
equation QBJ that can be cast into the form 

d t R{t) = M R(t) + I. (12) 

Here M. represents a generalized 15 x 15 Bloch matrix, 
the vector J" is an inhomogeneity with components 

I = (0, 0, 0, 0, 0, 0, 0, i 0, 0, 0, 0, 0, 0, -i IT, 0)* (13) 

and the vector R contains the matrix elements 
Qij = ( i \ q\ 3 ) °f t ne density operator, 

R = (Qll,Ql2,Ql3,QU,Q21,Q22,§23,Q2i, (14) 
031, 032, £?33, 034, 041, 042, 043) ■ 

The stationary solution of Eq. (|12(l is formally given by 
R st = -M- X I , (15) 
and an evaluation of the latter equation yields 
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The remaining non-zero components of R st are deter- 
mined by 

011 = 022 , 031 = 013 and 042 = 024 • (17) 
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In the case of the degenerate system, the population 
of the two ground levels will be equal and we have 
g>i3 = — £>24- Note that the minus sign arises since the 
dipole moments d\ and di are anti-parallel, and the co- 
herences ^14 and 023 are equal to zero because the driving 
field does not couple to the a transitions. 

In this paper we focus on the total intensity and the 
spectral distribution of the fluorescence light emitted by 
the atom in steady state. The total intensity 

I st = (E { -\r,t)-E i+ \r,t)) st (18) 

is given by the normally ordered first-order correlation 
function of the electric field, and the spectrum of reso- 
nance fluorescence is determined by the Fourier trans- 
form of the two-time correlation function of the electric 
field [H, 

OO 

S(Lu) = h J ^ (E { -\r,t + T)- E i+ \r,t)) st dT . 

— OO 

(19) 

In these equations, E [e + ^ denotes the negative 

(positive) frequency part of the electric field operator. 
At a point r = rf in the far-field zone, the negative 
frequency part of the electric field operator is found to 
be 

E { -\r,t) = E { J c (r,t) (20) 

4 

-^^^ 2 rx(rxd 4 )^ + (t)e^', 



where t — t — - is the retarded time, r\ = 1/ (AtteqC 2 ) 
and Sf — exp(=Fzwii) Sf. The first term stands for the 
negative frequency part of the free field. It does not con- 
tribute to the normally ordered correlation functions in 
Eqs. i|18[l and i|19[l as lon g as the point of observation lies 
outside the driving field 39] . The second term describes 
the retarded dipole field generated by the atom situated 
at the point of origin. 

Throughout this paper we assume that the point of ob- 
servation lies in the y-direction, where the z- and x-axes 
are defined by the polarization and the direction of prop- 
agation of the laser beam, respectively. An evaluation of 
the cross products in Eq. i|20|) shows then that the light 
emitted on the 7r transitions is linearly polarized along 
e z , whereas the light emitted on the a transitions is lin- 
early polarized along e x . The advantage of this detection 
scheme is that one can easily discriminate between the 
light emitted on the 7r and a transitions by means of a 
polarization filter. For this reason we will discuss the 
fluorescence light of the 7r-and a transitions separately. 



III. SPECTRUM OF RESONANCE 
FLUORESCENCE - tt TRANSITIONS 

We begin with a brief discussion of the steady-state in- 
tensity recorded by a broadband detector that observes 
the light emitted on the it transitions. According to 
Eqs. (HHJ and (HOJ), we have 

2 

il = <t>* E Hj{St~SjU, (21) 

where it was assumed that lo\ ~ LU2 to obtain a com- 
mon prefactor 0^ that we set equal to one in the 
following. The terms 7^ are defined in Eq. QlOp. 
and 712 = 721 = — \/7i72 describe the cross-damping be- 
tween the 7r transitions that arises as a consequence of 
quantum interference. However, these interference terms 
do not contribute to the total intensity, regardless of what 
the steady state solution might be, because the ground 
states are orthogonal, 

(£+5 2 -) st = (|l}(3||4)(2|) st = 0. (22) 

Consequently, the intensity emitted on the it transitions 
is not altered by interference terms and simply propor- 
tional to the population of the excited states, 

7- = ivy(en + 022) ■ (23) 

We now turn to the the spectrum of resonance fluores- 
cence emitted on the ir transitions. With the help of 
Eqs. (|T§j) and (|2U|) we arrive at 

2 00 

S*(aj) = \Y. Re / e ~ tGjT ^t(i + (*) >«t dr , 
i,j=l 

where G) = u> — lul is the difference between the ob- 
served frequency and the laser frequency. In contrast 
to Eq. 122|) . the terms proportional to 712 are now de- 
termined by the two-time averages (S^~(i+ T)S' 2 _ (t) ) s t 
rather than by the one-time averages. Indeed, we find 
that the correlation function 

Gi 2 (r)^-Vl^(S+(i + T)S^(i)) st (25) 

is different from zero for r > 0; a plot of Gn is shown 
in Fig. But this implies that there is quantum in- 
terference in the spectrum of the light emitted on the 7r 
transitions, although there is no interference in the to- 
tal intensity. To illustrate this result we decompose the 
transition operators in Eq. I)25|l in mean values and fluc- 
tuations according to 

St = (St) st l + 5S± . (26) 

The correlation function Gi2(r) becomes then 

Gi 2 (r) = -y/yw[(8St(i + T)5S^{i)U ( 27 ) 

+(5+) st (^)st]. 
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FIG. 2: Plot of the correlation function G12 in relation to 
its long-time limit Gi2(oo) = — - v /7i72( } Bt { §2 ) at for the 
degenerate system. The parameters are Q, — 3 x 10 s" 1 , 
A = 5 x 10 6 s _1 and 7 = 10 7 s _1 . G12 has to vanish at r = 
since the ground states are orthogonal. 



The two-time average of the fluctuations can be calcu- 
lated from the generalized Bloch equations and the quan- 
tum regression theorem (see Appendix) . It decays expo- 
nentially with a time constant on the order of 7 -1 and 
does not contribute to G\i in the long-time limit r — > 00. 
The mean values (Sf)^ — Q31 and (Sj ) s t — Q42 are 
given by matrix elements of the steady-state density- 
operator in Eq. (|l(jfl and are both different from zero. 
This is obvious from a physical point of view since the 
laser field creates a coherence on both transitions 1 <-> 3 
and 2 ^> 4. Consequently, the long-time limit of G\i 
reads 6*12(00) = — ^7172 ( Si )st( ) s t- It follows that 
the interference terms will affect the coherent and inco- 
herent spectrum of resonance fluorescence. 

Before we give expressions for the spectral distribution 
of the emitted light, we calculate the respective contribu- 
tions of coherent and incoherent scattering to the inten- 
sity To this end we apply the decomposition of the 
transition operators Eq. i|2rj[l to Eq. I|21|) . This allows us 
to write I£ t as the sum of four terms, 

-Ct = 4°oh + -Ch + ^inc + lint ■ ( 28 ) 

The first two terms account for the contribution of co- 
herent scattering (subscript "coh") and are given by 

4V = 7i|<Si + )st| 2 +7 2 |<4 + >st| 2 (29) 
41 = "2 VtT^ Re< S+ ) st ( S 2 " ) st . (30) 

In this equation, /° oh stands for the contribution of terms 
proportional to 711 and 722, and /™h i s the weight of 
the interference terms that can be positive or negative. 
By contrast, the sum of 7° oh and is the weight of 
the Rayleigh line that is always positive. The last two 
terms in Eq. (|28|l denote the contribution of incoherent 
scattering (subscript "inc"), 

IL = 7i ( SS+SS; ) st + 72 < SS+SS^ ) st (31) 
Il? c = -2^^Re(6S+SS^) st . (32) 



Since the ground states are orthogonal, Eq. I|26() allows 
to establish the relation 

( 6S+6S2 ) st = - ( S+ ) at < 5 2 - ) st . (33) 

If this expression is applied to Eq. (|32|) . it follows from 
Eq. H3Ufl that the interference terms J™^ and Jj"* are of 
opposite sign, i.e. 

41 = -tit ■ (34) 

This relation clarifies that the interference terms alter 
the weights of the coherent and the incoherent part of 
the spectrum, whereas the total intensity remains un- 
changed. Note that this is in contrast to the V-system 
with non-orthogonal transition dipole moments men- 
tioned in the introduction, where both the fluorescence 
spectrum and the total intensity show a signature of in- 
terference Hum 13. 

We now turn to the spectral distribution of the 
fluorescence light and employ Eq. (|2l)|l to write the 
spectrum of resonance fluorescence in Eq. Ij24(l as the 
sum of the coherent and the incoherent spectrum, 
S^Q)=S- oh (Q) + S^ c (uj), where 

S^(Q) = (/ c ° oh + 4oh)^) (35) 

2 00 

\ E ^ Re / ^ T (SS+(i+r)SSj(i)) st dr. (36) 

i,j=l o 

These two contributions will be discussed in the following 
Sections. 



A. Coherent spectrum of resonance fluorescence 

The coherent part of the fluorescence spectrum consists 
of the Rayleigh peak centered at u) = lul- In order to get 
a better understanding of how the weight of this line is 
affected by interference, we write it as 

4V + i l cl = I VTi( st ) 8t - Vts < s+ U I 2 • (37) 

In this equation, ( ) s t is proportional to the scattering 
amplitude on the 1^3 transition and — (§2 )st corre- 
sponds to the scattering amplitude on the 2^4 tran- 
sition. Note that the minus sign arises since the dipolcs 
di and d 2 are anti-parallel. Depending on the relative 
phase and the absolute values of the coherences ( S± ) s t 
and ( §2 )st, there will be constructive or destructive in- 
terference in the coherent part of the spectrum. We will 
now demonstrate that the degree of interference in the 
coherent spectrum can be controlled by means of the dif- 
ference 5 between the resonance frequencies of the it tran- 
sitions. Therefore, we write Eq. (|37(l as 

^coh + 4oL=^oh[l + C] , (38) 
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FIG. 3: Plot of the relative weight of the interference 
terms C(5) for different values of the detuning A of the 
laser field from the 1 <-> 3 transition. The parameters 
are given by 7 = 10 7 s _1 , A = —4 x 10 7 s _1 (solid line) and 
A = -5 x 10 6 s _1 (dashed line). 



where C — ^cSh/^coh ^ s the relative weight of the inter- 
ference terms. An explicit expression for C can be found 
with the help of the definitions in Eq. (|30B and the steady- 
state solution for g in Eq. Hlt)|) . 



C 



7 2 /4 + A(A-S) 
7 2 /4+,5 2 /4+(A-(5/2) 2 



(39) 



The absolute value of this quantity can be regarded as the 
degree of interference in the coherent spectrum. Figure|3| 
shows a plot of C as a function of S for two different 
(negative) detunings A. It is evident that C is equal to 
one in the case of the degenerate system. Therefore, we 
have perfect constructive interference for 5 — 0. In this 
case, the detunings on both n transitions are equal and 
hence we have (>5q ) s t = — (S2 )st, the two transitions 
are now perfectly equivalent. In addition, the weight 
of the Rayleigh line is then, apart from the branching 
probability b^, identical to the corresponding expression 
for a two-level atom |40| . 

As \5\ increases, C(S) decreases monotonously and be- 
comes zero at <5o = A[l + 7 2 /(4A 2 )] . Note that So can 
be either positive or negative, depending on the sign of 
A. In the case of A 2 ^ 7 2 , we have 60 ~ A. This im- 
plies that the interference term vanishes if the laser field 
is resonant with the 2 <-> 4 transition. The minimum 
of the curve is reached at <5 m i n = 2A(1 + 7 2 /(4A 2 )) and 
given by C(S min ) = -1/(1 + 7 2 /(2A 2 )). Consequently, 
C(^min) tends to —1 provided that A 2 ^S> 7 2 . The weight 
of the Rayleigh peak becomes then zero as a consequence 
of destructive interference, and the emitted radiation is 
solely incoherent. Note that this situation occurs if the 
detunings on the 1 <-> 3 and 2 <-> 4 transitions are ap- 
proximately equal and of opposite sign. In this case, 
the coherences (S^ ) s t and (S^~) s t cancel each other in 
Eq. QZJ- Finally C tends to zero as \S\ becomes much 
larger than |A| and 7. This is due to the fact that 
the interference term in Eq. H3Ufl consists of the prod- 
uct of ( §1 ) s t and ( §2 )st- If the detuning on one of the 



two 7T transitions becomes very large, 7™^ tends to zero, 
whereas J° oh remains different from zero. 



B. Incoherent spectrum of resonance fluorescence 

It is possible to evaluate the expression for S[ nc in 
Eq. (|36[) analytically, an outline of the calculation can 
be found in the Appendix. However, the general result 
is too bulky to present it here. We just mention that the 
spectrum does only depend on the difference S between 
the Zeeman splittings of the ground and excited states, 
but not on the parameter B (see Fig. In the case of 
the degenerate system, we find 



SUP) 



7 7 



2\tt\ 2 +Cj 2 



2 7 |fi| 



* 7T 7 2 /4 + A 2 + 2|ft| 2 \P(- 



(40) 



where P(z) is a cubic polynomial as a function of z that 
is defined as 

P(z) = i(z + 7 )[(2z + 7) 2 + 4A 2 ] +2(2z + 7 )W 2 . (41) 

Apart from the branching probability 6 W , this result is 
exactly the same as the incoherent spectrum of resonance 
fluorescence of a two-level atom [4(| . 

As soon as S becomes different from zero, the inco- 
herent spectrum differs considerably from the two-level 
spectrum. This is demonstrated in Fig. 0] a) which dis- 
plays S? nc for <5 = (dashed line) and 5 = -4 X 10 6 s _1 
(solid line). For 5 ^ 0, an additional central peak occurs 
whose width is much smaller than the decay rate 7. 
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FIG. 4: Incoherent spectrum of resonance fluorescence ac- 
cording to Eq. 136H . Plot (a) shows for the degener- 
ate system (dashed line) and for 8 — —4 x 10 6 s _1 (solid line), 
the other parameters are 7 = 10 7 s _1 , A = —4 x 10 7 s _1 and 
11 = 6 x 10 6 s _1 . In (b) and (c) the values of 5 are given by 
S = 80 and 8 — <5min, respectively, the other parameters are 
the same than in (a). Plot (d) shows the incoherent spectrum 



for the set of parameters A = - 
7 = 10 7 s" 1 and 8 = -8 x lfjV 



5 x lCTs -1 , n = 6 
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FIG. 5: Dressed state analog of the bare state system in 
Fig. The frequency of the laser field is labeled by u>l- 
For <5 7^ 0, the detuning of the laser field will be different on 
each of the n transitions. There are thus two effective Rabi 
frequencies Sli and Q2 involved. The splitting of the dressed 
states for fixed N is not to scale. 



Section IIII Al provides a detailed discussion of the 
weight of the interference term J™^ in the coherent spec- 
trum. These results can also be applied to the weight 
of the interference term in the inelastic spectrum by 
means of Eq. I|34|) . For example, it follows that the weight 
of the interference term in the inelastic spectrum van- 
ishes for S = Sq- This situation is shown in Figure Elb), 
where the width and the weight of the additional peak 
is larger than in a). For S — (5 m in and the parameters of 
Fig. 01 we know from Sec. IIII Al that the weight of the 
Rayleigh line is approximately zero. The corresponding 
incoherent spectrum is shown in Fig. c). Instead of the 
elastic delta-peak in the coherent spectrum we thus have 
a very narrow peak that occurs in the incoherent spec- 
trum. 

Finally, Fig.^d) shows S^ nc for a strong laser field. In 
this case, the weight of the interference terms is negligible 
as can be verified with the help of Eq. (|30() . However, 
the incoherent spectrum still deviates from the Mollow 
spectrum if 5 =/= 0. This can be easily understood with 
the aid of the dressed states [38l l4l| of the system. If 
N denotes the number of laser photons of frequency lul, 
the dressed states can be expressed in terms of the bare 
states as follows, 

\1{N)) = e 10 sin 61 1 1, AT )+ cos 61 1 3, AT + 1 ) 

|2(A0) = e^cosOi 1 1, AT) — sin 61 1 3, N+l) (42) 

where tan 26 1 = 2 |f2|/A and 

|3(AT)) = e^sin© 2 |2, N) -cos6 2 |4, N+l) 

|4(A0) = cos 2 1 2, AT ) + sin 6 2 1 4, N+l) (43) 

with tan26 2 = 2 |0|/(A - S) (0 < 1; 9 2 < tt/2, e** = 
f2/|n|). Figure shows the relative position of the 



dressed states for two manifolds with N and (N—l) laser 
photons, respectively. Note that 1 1(AT) ) and | 2(AT) ) are 
separated by a frequency interval of Cli = ^/ 4| S~J | 2 + A 2 , 
whereas the spacing between 1 3(AT) ) and 1 4(AT) ) is 
given by f2 2 = v /4|0| 2 + (A - S) 2 . The sidebands in 
the spectrum of the ir transitions result from the tran- 
sitions |l(Af)>^|2(Ar-l)>, | 2(JV) > — 1 1(JV - 1) >, 
I 3(N) ) -> I 4(AT - 1) ) and | 4(A r ) > — > | 3(7V - 1) ). Con- 
sequently, they will be located at the frequencies u>l ± ^1 
and Wx, ± Q 2 . For 6 ^ 0, we thus expect four sideband 
peaks symmetrically placed around the laser frequency 
iv l, precisely as in Fig. H|d). 



C. Influence of the interference terms on the 
fluorescence spectrum 

In this Section we investigate how the interference 
terms alter the fluorescence spectrum emitted on the it 
transitions. Here we only consider the degenerate system 
that is distinguished by maximal constructive (destruc- 
tive) interference in the coherent (incoherent) part of the 
fluorescence spectrum, see Sec. Mil If the interference 
terms in Eq. (|2"4"|l are omitted, the fluorescence spectrum 
reads 



1 2 

S&(fl) = ~y>iHe 

7T e — ' 

i=l 



-(S+(t + T)S7(t)) st dr 



(44) 

The fluorescence spectra with and without the interfer- 
ence terms according to Eqs. I|24|) and (|44l) are shown in 
Fig. for different parameters of the laser field. If the 
saturation parameter defined in Eq. (|49(l is much larger 
than unity, the weight of the interference terms goes to 
zero. However, Fig. a) demonstrates that the interfer- 
ence terms still alter the shape of the fluorescence spec- 
trum in the region of the sideband peaks. The spectrum 
S 1 *' with interference terms is identical to the fluorescence 
spectrum of a two- level atom (see Sec Mlf) . and thus the 
ratio between the central and the sideband peaks reads 
1:3:1. For the spectrum without the interference 
terms and a branching probability of b n = 1/3, this ratio 
reads 7 : 15 : 7. 

FigureHJlb) shows S n and Sq for low saturation. In this 
case, the spectrum without interference terms is distin- 
guished by a narrow peak centered at the laser frequency 
that occurs in addition to the clastic Rayleigh peak. A 
numerical analysis shows that Sq can be written as 



sjp) « iUm + s? nc (u>) + s; eak (Q) . 



(45) 



In this equation, the first term represents the Rayleigh 
peak whose weight misses the interference term 7™^ that 
is present in Eq. I|35|l . The second term stands for the 
incoherent spectrum according to Eq. (|40|l . The last term 
describes a Lorentzian of weight J™^ and width T w that 
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FIG. 6: Fluorescence spectrum for the degenerate system 
according to Eq. 12411 . The solid line (dashed line) shows the 
spectrum with (without) the interference terms proportional 
to 712, 721- The Rayleigh peak (the vertical line at u — lul) is 
present both with and without interference terms. Note that 
its weight is larger if the interference terms are taken into 
account. However, the sums of the integrated coherent and 
incoherent spectra with and without the interference terms 
are identical, making the total intensity independent of the 
interference terms. In a), the parameters are Q. = 5 x 10 7 s -1 , 
A = and 7 = 10 7 s' 1 . For b), we have Q = 10 7 s"\ A = 
2 x 10 7 s _1 and 7 = 10 7 s _1 . 



is centered at the laser frequency, 



pnt 

i coh 

7T w 2 + r 2 



(46) 



The weight of the extra peak S^ eak is determined by the 
constraint that the total intensity is independent of the 
interference terms (see Sec. Illlf) . Therefore, Sp 0ak has to 
compensate for the reduced weight of the Rayleigh line of 
Sq as compared to the spectrum with interference terms. 
In general, the width r„. of the extra peak Sp Cak is smaller 
than the decay rate 7. If the saturation parameter s is 
much smaller than unity, we find (b n = 1/3) 



rint 
coh 



2s)s and r„ 



22(3- 



5s)s . (47) 



Figure EJb) allows to summarize the effect of the inter- 
ference terms on the fluorescence spectrum in the case 
of low saturation as follows. The spectrum without in- 
terference terms displays a narrow peak 5p Cak of finite 
width at the laser frequency that is absent if the interfer- 
ence terms are taken into account. Therefore, quantum 
interference cancels the incoherent response of the atom 
at the laser frequency lol- 

In conclusion, the experimental observation of the flu- 
orescence spectrum confirming the solid lines in Fig. 
would give evidence for vacuum-mediated interference ef- 
fects as described by terms proportional to 712. So far, 
interference effects of this kind have not been observed 
in atomic systems. 



IV. SPECTRUM OF RESONANCE 
FLUORESCENCE II - a TRANSITIONS 

This Section is concerned with a brief discussion of 
the fluorescence spectrum emitted on the a transitions. 
Since the laser field does not couple to these transitions, 
the spectrum contains only an incoherent part. We arrive 
at 



4 00 



where 



(t+T)SSr(i)) st dT , 
(48) 

is a geometrical factor that we set equal to one 
in the following. It has been pointed out in Sec. ITU that 
the light emitted on the a transitions is linearly polarized 
along e x if the point of observation lies in the y-direction. 
Therefore, the cross terms ( SS^(t + T)SS^(t) ) st and 
( 5S^ (i + r)SS^(t) ) st will, in principle, contribute to 
the spectrum in Eq. 148fl . However, we find that the 
latter two-time averages are equal to zero. For different 
driving schemes where the laser field couples to the a 
transitions, the cross-correlation terms have to be taken 
into account as is the case in the work of Polder et. 
al. j3^| . The exact analytical expression for S a is too 
bulky to display it here. Instead we will discuss S" in 
the case of the degenerate system (B = 5 = 0) and for 
different regimes of the driving field strength that will be 
characterized by means of the saturation parameter 



2|fi| 



A 2 + 7 2 /4 



(49) 



In the range from a weak to a moderately strong laser 
field (s < 1), a numerical analysis reveals that S a can be 
written as 



(50) 



In this equation, the first term stands for the incoher- 
ent spectrum of a two-level atom according to Eq. 140fl . 
The prefactor b a /b 7r accounts for the different branch- 
ing probability of the a transitions as compared to the tt 
transitions. The second term represents a narrow peak 
that is centered at the laser frequency u> = lol- It can be 
modeled as a Lorentzian of weight W a and width Y a , 



tt w 2 + r 2 



(51) 



The weight of Sp Cak is determined by the total intensity 
emitted on the er transitions 



baiim + Q22) 

and the weight of b^/b^ S^ c . We arrive at 
W a =4& CT7 |p 13 | 2 , 



(52) 



(53) 



where Qi 3 is given in Eq. ljTo|> . The width T a of the 
additional peak is smaller than the decay rate 7. If s 
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FIG. 7: Spectrum of resonance fluorescence emitted on the 
a transitions (solid line) in comparison with the fluorescence 
spectrum of a two- level atom (dotted line). The parameters 
in (a) are fi = 5 x 10 6 s~ x , A = 6 x 10 6 s" 1 and 7 = 10 7 s _1 . 
In (b), the parameters are f2 = 6 x 10 7 s -1 , A = and 
7 = 10 7 s _1 . Note that S a deviates slightly from the Mol- 
low triplet in the region of the sideband peaks. 



is much smaller than unity, the width and the weight of 
Speak are S iven by 



b a j{l-2s)s 



(54) 



At the same time, the contribution of S[ nc b^jb^ to S a is 
small such that the spectrum is dominated by the central 
narrow peak 5p eak . If the field strength is increased, the 
weight of the extra peak Sp Cak gets smaller. Figure a) 
shows S a (solid line) and S? nc b^/b^ (dotted line) for a 
moderately strong laser field, the saturation parameter 
is on the order of unity. Nevertheless, the spectrum S a 
is still dominated by the sharp peak S* eek that exceeds 
the central peak of the two-level spectrum by one order 
of magnitude. 

For a strong driving field (s> 1), the weight of S° cak 
goes to zero and the central peak of S° coincides with the 
corresponding peak of the Mollow spectrum. However, 
the sideband peaks of S" 7 differ from those of a two-level 
atom as can be seen from Fig. 0b). In the secular limit, it 
is advantageous to employ the dressed state picture in or- 
der to obtain analytic expressions for the sideband peaks, 
being well separated from the central peak whose analytic 
form can be taken over from the well-known results for 
a two-level atom [3^, ^(J . The fluorescence spectrum for 
a resonant driving field can be achieved by a tedious but 
straightforward calculation that follows the procedure of 
Chapter VI.E in [H, 



S°(Co) 



8tt 



(55) 



+ 7 



7 /2 



47T 7 2 /4 + o> 2 



+ 7 



• sb 



8tt F s 2 b + (fii + Co) 2 1 



where fi x = ^/4|f2| 2 + A 2 . A comparison of the latter 
equation with the corresponding expression for the Mol- 
low spectrum reveals that the weights of the sideband 



peaks differ only by the branching probability b„. For 
the width of the sideband peaks in Eq. I|55|) we find 



| = ~(3-& CT )7- (56) 



r S b = ^V7i72 



Note that the second equality is obtained by virtue of 
Eq. . The ratio between the heights of the central 
peak at Co = and the sideband peaks at Co = ±Qi is 
found to be 3 — b a . For b a = 2/3, the peak ratio is thus 
3:7:3. By contrast, the peak ratio of the Mollow 
spectrum reads 1:3:1. A precise measurement of the 
peak ratio would thus provide a means of determining 
the branching probability b a of the degenerate system 
experimentally. 

Note that the width of the sideband peaks in Eq. H5ti|) 
depends on the cross-damping terms ^pyCyi that appear 
in the master equation through the spontaneous emission 
term C 1 g in Eq. JSJl. If these interference terms were not 
present, the peak ratio would not depend on the branch- 
ing probabilities and would be given by 1 : 2 : 1. The 
spectrum emitted on the a transitions shows thus an in- 
direct signature of interference. 



V. DISCUSSION 

In Section |ml we have shown that the interference 
terms proportional to 712 contribute only to the spectrum 
of resonance fluorescence, but not to the total intensity 
in Eq. 12111 . Here we demonstrate that this result is a 
consequence of the principle of complementarity, applied 
to time and energy. 

If the total intensity is measured, complementarity 
does not impose any restrictions on the time resolution 
of the measurement since the photon energies are not ob- 
served. It is thus possible to observe the temporal aspect 
of the radiative cascade, i.e. one could determine the 
photon emission times. The time evolution of the driven 
atom is then most suitably described in the bare state 
basis. For example, assume that the atom is initially in 
ground state | 3 ) . The laser field will induce Rabi oscilla- 
tions between the excited state | 1 ) and | 3 ) . Immediately 
after the spontaneous emission of a photon, the atom is 
found in ground state | 3 ) (it transition) or | 4 ) (a tran- 
sition). Subsequently, this sequence of Rabi oscillations 
and a spontaneous emission event is repeated. In this 
description, each emission process on one of the it transi- 
tions is independent of the other it transition. However, 
quantum interference does only occur if various indistin- 
guishable transition amplitudes connect a common initial 
state to a common final state. Since the tt transitions do 
neither share a common initial nor a common final state, 
we must conjecture that the total intensity is not affected 
by interference. 

The lack of interference in the total intensity can also 
be explained by drawing an analogy to the two-slit ex- 
periment. It is well known that the interference pat- 
tern vanishes as soon as it is principally possible to know 
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through which of the two slits each object (electrons or 
photons) has moved. Similarly, the internal states of our 
atom can be regarded as a which- way marker. Since the 
experimental conditions allow, at least in principle, to 
determine the atomic ground state immediately after the 
detection of a 7r-photon, one could decide on which of the 
two 7r-transitions the photon was emitted. Consequently, 
the observer could reveal the quantum path taken by the 
system and hence there is no signature of interference. 
Note that this argument requires that the retardation be- 
tween the times of emission and detection is much smaller 
than the time between successive emissions. This condi- 
tion can typically be achieved in atomic systems. 

A totally different situation arises if the detector mea- 
sures the spectrum of resonance fluorescence and hence 
the energy of the emitted photons. First of all, it is ad- 
vantageous to illustrate the energy aspect of the cascade 
of spontaneously emitted photons in the dressed state 
picture (see Fig. |HJ rather than in the bare state pic- 
ture. The crucial difference between the measurement of 
the total intensity and the fluorescence spectrum is the 
following. In the latter case, the observer decides to de- 
termine the photon energies precisely. Since time and en- 
ergy are complementary observables, the temporal aspect 
of the radiative cascade is not accessible simultaneously. 
Next we demonstrate that precisely this lack of informa- 
tion about the temporal sequence of photon emissions 
allows for the interference mechanism in the fluorescence 
spectrum. 

A quantitative description of time-energy complemen- 
tarity is achieved via the time-energy uncertainty. If the 
photon energies are determined with a precision of Acj, 
the time-energy uncertainty relation enforces that the 
time of observation has to be at least on the order of 
1/Aw. Since the observer can only notice the detection 
of a photon after the observation time has elapsed, the 
photon emission times are indeterminate within a time 
interval of At = 1/Aw. For the moment we envisage an 
ideal measurement of the fluorescence spectrum. In this 
case, the atom will emit (infinitely) many photons dur- 
ing the (infinite) time of observation. In addition, the 
photon emission times are indeterminate, and thus the 
time order in which these photons have been emitted is 
unknown. It follows that the transition amplitudes cor- 
responding to the various time orderings of the photons 
will interfere. 

We illustrate this interference mechanism on the basis 
of Fig. that shows a cascade of only two photons, one 
emitted on a tt transition and the other on a a transi- 
tion. Assume that the atom is initially in the dressed 
state | 4(7V) ). In one of the two cascades, the atom de- 
cays first to the state | 4(iV — 1) ) by the emission of a 
tt photon on the bare state transition | 2 ) — > | 4 ) . The 
subsequent emission of a er photon takes the atom to the 
state 1 1(N — 2) ) within the manifold with N — 2 laser 
photons. In the second cascade, the time order of the 
two photons is reversed. The atom decays now first to 
the state 1 1(N — 1) ) by the emission of a a photon, and 
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FIG. 8: Radiative cascade in the dressed states of the system 
[see Eqs. 142H and 1431 1. The splitting of the dressed states 
for a fixed number of laser photons N is not to scale. Each of 
the two indicated cascades involves the emission of a n pho- 
ton and a a photon with wave vector fe^ and k a , respectively. 
Depending on the time order of their emission, the tt pho- 
ton is either emitted on transition | 4(7V) ) — > | 4(7V — 1) } or 
| 1(N — 1) } — > | l(iV — 2) }, corresponding to the bare state 
transitions 2 } — > 4 } and 1 } — > | 3 } , respectively. Since the 
final and initial states of the two cascades are identical, the 
corresponding transition amplitudes may interfere. 



then to the final state 1 1(N— 2) ) under the emission of a 
tt photon. In contrast to the first cascade, this tt photon 
is now emitted on the bare state transition 1 1 ) — > | 3 ) . 
Since the two cascades in Fig. [5] have the same initial 
and final states, and since it is in principle impossible 
to determine the quantum path taken by the system, the 
two transition amplitudes corresponding to different time 
orders of photon emissions interfere. In one of the transi- 
tion amplitudes the tt photon stems from the | 2 ) — > | 4 ) 
transition, and in the other from the 1 1 ) — * | 3 ) transi- 
tion. Exactly this mechanism gives rise to the interfer- 
ence effects in the fluorescence spectrum that are medi- 
ated by the cross-damping terms in Eq. H24JI . 

The provided explanation can also be employed to il- 
lustrate why there is no interference in the fictitious sit- 
uation of perpendicular dipole moments d\ and d^. In 
this case, a photon can either stem from d\ or d2, but 
not from both transitions. It is then impossible to realize 
both cascades in Fig. |S| and hence there is no interfer- 
ence. Moreover, it becomes now clear why the spectrum 
emitted on the a transitions depends on the interference 
terms 712 and 721. For anti-parallel dipole moments d\ 
and c?2 there are two transition amplitudes that involve 
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the emission of a a photon , and for perpendicular dipole 
moments there would be only one. We emphasize that 
the discussion has been restricted to a cascade of only 
two photons for the sake of simplicity. In principle, all 
possible cascades with an arbitrary number of photons 
have to be considered, but the general idea remains the 
same. 

It is also possible to provide an explanation for the in- 
terference in the coherent spectrum, but the elastic scat- 
tering events cannot be visualized in the dressed state 
basis. However, in the case of low saturation (s <C 1) 
the process of elastic scattering can be illustrated in the 
bare state basis such that the atom hops from one ground 
state to another by the absorption of a laser photon and 
the emission of a scattered photon. The excited states 
act as intermediate states and can be adiabatically elim- 
inated. Since it is impossible to tell on which of the two it 
transitions the photon was scattered, it is plausible that 
one has to sum the scattering amplitudes first and then 
take the absolute value squared in order to obtain the 
weight of the Raylcigh line in Eq. (|38H . 

Next we demonstrate how the interfering transition 
amplitudes that correspond to different time orders of 
photon emissions enter the expression for the spectrum 
of resonance fluorescence in Eq. Q24[). Let (aj.) be the 
annihilation (creation) operator of a photon in a mode 
of the radiation field that is actually observed by the de- 
tector, being sensitive only to photons emitted on the 7r 
transitions. The rate at which the photon number in this 
particular mode changes is given by 



R n (t) = d t (at(t)a 7r (t)] 



(57) 



If one follows the lines of Chapter 7 in 2] , one can show 
that the steady-state value of is proportional to the 
spectrum of resonance fluorescence, 



lim R v {t) ~ S*(c\K\ 



(58) 



In this equation, k n denotes the wave vector that corre- 
sponds to the observed mode a w , and c is the speed of 
light. In order to evaluate the left hand side of Eq. 1)58(1 . 
we will label the basis states \i(N); {n}) of the total 
system (atom + laser field + vacuum modes) by three 
quantum numbers, namely the dressed states i, the num- 
ber of laser photons N and the state of the vacuum modes 
{n}. The mean value on the right hand side of Eq. I|57|) 
becomes then 



4 

E 

= 1 N, {n} 



where \C l N ( n \(t)| 2 is the probability to find the system 
at time t in state | i{N); {n} ) and N n ({n}) is the expec- 
tation value of a\a-K in this state. We assume that the 
system is in some initial state | V'o ) at time t = with all 
vacuum modes being empty. If the time evolution oper- 
ator is labeled by U(t, 0), the transition amplitude from 



the initial state | ipo ) to the final state \ i(N); {n}) can 
be written as 



C l N , {n} (t) = (i(N); {n}\U(t,0)\^o] 



(60) 



Let us assume that the final state contains q scattered 
photons that are characterized by their wave and po- 
larization vectors, {n} — {k^e^, ^2^2, • ■ ■ , k q e q }. We do 
not attempt to evaluate Eq. (|(jl)fl explicitly, but in prin- 
ciple one would introduce q — 1 intermediate states and 
arrange the q scattered photons into a certain order. But 
since there is no distinguished time order of the scattered 
photons, there are, in principle, q\ transition amplitudes 
involved in the evaluation of Eq. H6U|) that will all inter- 
fere. 

In conclusion, we demonstrated that the interference 
in the spectrum from the 7r transitions can be explained 
in terms of interference between transition amplitudes 
that correspond to different time orders of photon emis- 
sions. If the spectrum of resonance fluorescence is ob- 
served, the principle of complementarity enforces that 
these transition amplitudes are indistinguishable. If the 
total intensity is recorded by a broadband detector, the 
temporal aspect of the radiative cascade can in principle 
be observed. Consequently, the possibility of interfer- 
ence between different time orders of photon emissions 
is ruled out. The preceding discussion of our results 
also implies that the experimental setup — potentially af- 
ter the photon emissions — decides if interference takes 
place, a feature that is also known from quantum eraser 
schemes pHH^ . 

We now refine our analysis and consider a detector 
with a finite frequency resolution Auj that allows us to 
study the continuous transition from perfect frequency 
resolution to perfect time resolution. For simplicity, we 
consider only the degenerate system (B = 5 = 0). If a 
filter of bandwidth A and setting frequency to is placed 
in front of a broadband detector, the spectrum can be 
determined with an accuracy of A, and the temporal res- 
olution is on the order of A -1 . The spectrum of resonance 
fluorescence emitted on the 7r transitions reads then [SI 



S*(w,A) = 



(61) 



2 



(st(i+r)s-(i)) st dr 



In the absence of interference terms the spectrum will be 
denoted by Sq(uj,X) and is obtained from Eq. (|61|l by 
omitting the terms proportional to 712 and 721. For the 
rest of this Section we assume that the saturation pa- 
rameter s is much smaller than unity. To a first approxi- 
mation, the incoherent contribution to the spectrum with 
interference terms can then be neglected. In the presence 
of the filter, the coherent 5-peak becomes a Lorentzian 
of width A and weight 7° 



^coh' an d t nus we obtain 
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FIG. 9: The solid lines show the fluorescence spectra recorded 
with a finite frequency resolution A. The dashed curves are 
the spectra without the interference terms proportional to 
712, 721 in Eq. 16111 . The parameters are f2 = 7 x 10 6 s _1 , 
A = 2 x 10 7 s"\ 7 = 10 7 s _1 and B = 8 = 0. This corresponds 
to a saturation parameter of s = 0.235 and a mean number 
of photons per unit time of approximately 9.4 x 10 5 s _1 . The 



filter bandwidths are given by a) A 
c) A = 1.9 x 10 6 s _1 and d) A = 10 7 



1(T 

-i 



-\ b) A= 10 4 



Similarly, we neglect the contribution of S[ nc to the spec- 
trum without interference terms in Eq. (|45|l . the (5-peak 
becomes a Lorentzian of width A and weight I® oh , and 
S'p Cak is replaced by a Lorentzian of width + A and 

weight 1™, 



SJ(<D,A)«i-sk- 
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rint 
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+ (IV + A) 2 



(63) 



Figure i shows the fluorescence spectrum according to 
Eq. I|61|l (solid lines) for different values of the filter band- 
width A and for low saturation. The dashed lines are 
the spectra without the interference terms. In Fig. a), 
the bandwidth A is much smaller than T„. Therefore, 
the widths of the lines S*(Cj, A) and Sq(Cj, A) are clearly 
distinct. If A is increased, the differences between the 
spectra with and without the interference terms dimin- 
ish until both curves are virtually identical for A = 7 
(Fig. id)). 

These results can be understood as follows. With an 
increasing filter bandwidth A, the smallest time interval 
At that can be resolved by the detector without violating 
the time-energy uncertainty gets shorter. Therefore, the 
observer can in principle obtain more information about 
the quantum path taken by the atom. Consequently, 
we expect that the signature of interference in the fluo- 
rescence spectrum diminishes for increasing A. This is in 
agreement with Fig.^and completely analogous to a two- 
slit experiment, where the visibility of the interference 
pattern is reduced at the cost of which-path information 
and vice versa . 

Furthermore, we demonstrate that the time-energy un- 
certainty relation allows to estimate the smallest filter 



bandwidth A for which the spectra with and without in- 
terference terms should be indistinguishable. Since the 
total number of photons emitted per unit time is equal 
to 7(^11 + ^22), the mean time between successive photon 
emissions is determined by = l/[7(gn + £"22)]- If the 
bandwidth A is chosen such that the temporal resolution 
could be much better than the mean time between succes- 
sive photon emissions, i.e. A -1 «0 = 1/ [7(011 + Q22)], 
we have 



A > 7(f?n + £22) ~ (1 - s)s^/2 



(64) 



Under these conditions, the radiative cascade of photons 
could be observed in a time resolved way and it is ex- 
tremely unlikely that more than one spontaneous emis- 
sion takes place during the time of observation. Since 
this rules out the interference mechanism as described in 
Sec. [V] the signature of interference in the fluorescence 
spectrum should disappear. But if inequality l|64|l holds, 
it follows that A ^> IV, and in this case S*(u>, A) and 
Sq(u>, A) are indeed indistinguishable as can be seen from 
Eqs. (S2 )l and This is confirmed by Fig. id) that 

shows Sq(oj, A) and S v (Cj,X) for a bandwidth A that is 
about ten times larger than the mean number of photons 
emitted per unit time. The two spectra are now virtually 
indistinguishable. 

It remains to explain the sharp peaks in the incoherent 
spectrum. To this end we return to Fig. 0] that shows 
the incoherent spectrum S[ nc for several values of the 
parameter S. A narrow central peak occurs only in case 
of the non-degenerate system (S ^ 0), and thus only if 
the weight of the Rayleigh line deviates from its maximal 
value attained at S = 0. Therefore, the narrow central 
peak in the incoherent spectrum may be regarded as a 
(partially) broadened coherent peak. This broadening 
can be understood as follows. Except for 5 = 0, the 
two 7T transitions are not equivalent since the absolute 
value and the phase of the coherences ( ) st and ( ) st 
will be different. The time that the atom spends on the 
1 «-> 3 transition can thus be regarded as a dark period 
with respect to the 2 <-> 4 transition and vice versa. This 
suggests that the sharp peaks in the incohere nt sp ectrum 
can be explained in terms of electron shelving |4a , l47ll4q . 
This explanation is also applicable to the sharp peak in 




FIG. 10: Schematic representation of elastic scattering events 
into the 3 — > 1 — > 4 (solid arrows) and 4 — » 2 — » 3 channels 
(dotted arrows). These processes account for the sharp peak 
in the fluorescence spectrum S" 7 emitted on the a transitions. 
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the spectrum from the a transitions. FigurelTUIillustrates 
the scattering events that give rise to this peak. If the 
atom is initially in state | 3 ) , a scattering event can bring 
it to state | 4) (solid arrows). The scattered photon has 
then been emitted on one of the a transitions. Before the 
next photon can be scattered on that same transition, 
a similar scattering process has to take place into the 
4^2-^3 channel (dotted arrows) . Consequently, every 
emission on one of the a transitions is followed by a dark 
period on that same transition. 

It should be mentioned that related interference effects 
between transition amplitudes corresponding to different 
time orders of photon emissions do also play a role in the 
fluorescence spectrum of other systems |49j. However, 
the distinguished feature of the system presented here is 
that this mechanism gives rise to interference effects be- 
tween the two 7r transitions that do not share a common 
state. 

We also point out that our system belongs to a class 
of setups that display interference and complementarity 
in the time-energy domain. In a conventional double-slit 
experiment, spatially separated pathways result in an in- 
terference pattern in position space. This is in contrast 
to our setup, where different temporal paths lead to in- 
terference in the energy domain. The work presented 
here is thus related to recent double-slit experiments in 
the time-energy domain [HcL l5lj | . In these experiments, 
ultra-short laser pulses of atto- or femtosecond duration 
open different time windows for the photoionization of 
an atom. If the energy spectrum of the photoelectrons 
is measured, these time-slits are indistinguishable and an 
interference pattern is observed. Moreover, it has been 
demonstrated that interference in the time-energy do- 
main can occur in the intensity correlations of different 
spectral components in a two-level atom jFJ H2, 13 • 



VI. SUMMARY 

The key result of this paper is that there is quantum 
interference in the spectrum of resonance fluorescence un- 
der conditions of no interference in the total intensity, 
being enforced by the principle of complementarity. For 
the system considered here, it claims that it is impossi- 
ble to observe the temporal and the energy aspect of the 
radiative cascade of the atom at the same time. If the 
fluorescence spectrum is observed, the photon emission 
times are indeterminate. The interference in the fluores- 
cence spectrum can thus be explained in terms of inter- 
ferences between transition amplitudes that correspond 
to different time orders of photon emissions. 

It has been shown that the degree of interference in the 
fluorescence spectrum emitted on the 7r transitions can 
be controlled by means of an external magnetic field. In 
particular, the degree of interference in the coherent part 
of the spectrum can be adjusted from perfect construc- 
tive to perfect destructive interference. Under condi- 
tions of perfect destructive interference, the weight of the 



Rayleigh line is completely suppressed. If the difference 
8 between the resonance frequencies of the 7r transitions 
is different from zero, the incoherent spectrum emitted 
on the 7r transitions contains a very narrow peak whose 
width is smaller than the decay rate 7. This peak has 
been identified as a partially broadened coherent peak 
and can be explained in terms of electron shelving. 

The spectrum emitted on the a transitions contains 
only an incoherent part. In the case of a weak driv- 
ing field and for the degenerate system, the fluorescence 
spectrum displays a narrow peak that can be regarded 
as broadened coherent peak. For a strong driving field, 
the widths of the sideband peaks differ from the Mollow 
spectrum. We have shown that the ratio between the 
peak heights of the central and the sideband peaks dis- 
play an indirect signature of interference. In addition, a 
measurement of the relative peak heights allows to de- 
termine the branching probability b a of the spontaneous 
decay of each excited state into the a channel. 
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APPENDIX: CALCULATION OF THE 
TWO-TIME AVERAGES 

In this section we outline how the functions 

00 

Sij(w) = Re J e-^ T {8S+ {t + t)5SJ {t) ) st dr . (A.l) 


can be evaluated by means of the quantum regression 
theorem [5^, |5|J . To this end we introduce the operators 
Aij that are connected to the atomic transition operators 
Aij (taken in the Schrodinger picture) by 

Aij = W^AijW , (A.2) 

where the unitary transformation W is defined in Eq. J5J • 
In particular, the operators & introduced in Sec. (|TT|l can 
be identified with the operators Aij according to 

S+ = A 13 S+ = A 24 S+=A 23 S+ = Au. (A.3) 

The corresponding Heisenberg operators are then defined 

as 

Aij{t) = U\t,0)AijU(t,0) , (A.4) 

and the time evolution operator has been labeled by U. 
A straightforward calculation shows that the mean values 
of the these Heisenberg operators are directly related to 
the matrix elements of the reduced density operator g in 
the rotating frame, 

(Aij(t)) =Tr A [AijQ(t)} = Qji(t) ■ (A.5) 
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In this equation, TrA [ ■ ] denotes the trace over atomic 
degrees of freedom. Next we arrange the operators Aij 
in a column vector 

L = (^11,^21, -431,^41,^12, ^22,-^32,^442, 

•4 13 , ^23 , -^33 , -^43 , -4l4 , A24 , -^34 ) * 

such that (L(t)) coincides with the Bloch vector R(t) of 
Eq. (JUJl, i.e. (L(t)) = R{t). It follows that the mean 
values (L(t)) obey the generalized Bloch Equation 1)120. 
If we decompose each component of L in mean values an 
fluctuations according to Aij = 6 Aij + (Aij ) s tl, we can 
cast ( L ) into the form 

(L(t)) = (6L(t))+(L) st , (A.6) 

where (L) st = R s t = -M~ l I. If Eq. (|A.6jl is plugged 
into Eq. JEJ we obtain a homogeneous equation of mo- 
tion for the fluctuations, 

d t (SL(t)) =M{SL(t)) . (A.7) 

The two-time correlation functions ( 5Li(t+T)6Lj(i) ) for 
i £ {1, . . . , 15} and fixed j can be written in vector no- 
tation as gi(t, t) — (5L(t + T)8Lj(t)). According to the 
quantum regression theorem, g 3 obeys the same equation 
of motion than the one-time averages ( 5L(t) ), 

d T g 3 = M g j for r > . (A.8) 



If G 3 (t, z) denotes the Laplace transform of g 3 (i, r) with 
respect to r, it follows 

G j (t,z) = [zl-M]~V(i,0) . (A.9) 

We need the Laplace transform at z — iuj in steady state 
to determine the functions Sij(uj) of Eq. (|A.1() . With the 
definitions 

R 3 = Jim g 3 (i, 0) and K 3 {Cj) = Jim G 3 (i, z = iw) 

t— >OC t — >oo 

(A.10) 

we arrive at 

K 3 (Q)= [iu>i- M^R 3 . (A.ll) 

The relevant correlation functions that are needed for the 
evaluation of Eq. (|86(l and l|48(l are then given by 

Sii(w) = Re [K 3 (w)] g S 21 (Q) = Re [K 3 (Q)] 14 

S 22 (Cj) = Re [K 8 (u)] M S 12 (w) = Re [K 8 (Cu)] g 

S 33 (Cj)=Re[K 7 (u)] 10 5 44 (^) = Re[i<: 4 (^)] 13 . 

(A.12) 

Finally, we remark that Eq. IjtjlH can be evaluated if one 
replaces iui in Eq. (| A. 1 1|> by iu) + A. 
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